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^ '. Abstract 
oo 

An accurate description of nuclear matter starting from free-space nuclear 

^ ' forces has been an elusive goal. The complexity of the system makes approx- 

22 ■ imations inevitable, so the challenge is to find a consistent truncation scheme 

^^ . with controlled errors. Nonperturbative effective field theories could be well 

^^ , suited for the task. Perturbative matching in a model calculation is used 

Q>^ I to explore some of the issues encountered in extending effective field theory 

0^ ' techniques to many-body calculations. 






I. INTRODUCTION 

Nuclear forces have been studied in depth over the past fifty years, leading to excellent 
phenomenological descriptions |l[]. A main goal of this enterprise has been to arrive at 
the "best" two-nucleon potential and to use it in many-body calculations to determine the 
properties of nuclear matter and finite nuclei. Effective field theory (EFT) analyses offer 
a different perspective: There is no best two-body potential. Indeed, the off-shell behavior of 
a potential or amplitude is not observable and should not influence the predictability of the 
theory if a consistent power-counting scheme is used 0. In addition, many-body contact 
terms are an inevitable consequence of using low-energy degrees of freedom and so must be 
included along with two-body interactions to eliminate off-shell ambiguities 0,§]. Finally, 
the symmetries of the underlying theory of QCD can be used to constrain the interactions, 
such as using chiral symmetry to describe two-pion exchange [^,0. These features merit a 
reanalysis of nuclear matter within the context of an EFT. 

Effective fleld theory techniques may also offer special advantages for many-body sys- 
tems, where exact solutions are not realizable. Estimates of the errors in approximate 
calculations and clear guidelines for making systematic improvements are hallmarks of an 
EFT. Even with a well-behaved underlying theory, like QED, an EFT can be beneflcial for 
organizing the contributions from various effects, as NRQED does for the radiative correc- 
tions to positronium P| . Similar beneflts may result from an application of EFT to nuclear 
matter. 

Skyrme wrote down the first effective theory for nuclear matter and used properties 
of nuclei to fit parameters in a mean-field lagrangian. Early attempts to motivate these 
parameters from the free-space nucleon-nucleon potential were promising [10], but were 



shown to depend on the form of the two-body potential [|11|. More recent studies have 
focused on scaling properties near the Fermi surface |T^-p!^ to determine which operators 
are important for Fermi liquids |jl5| as well as to describe the pairing gap [|l^. The behavior 



of nuclear binding in the chiral limit |T^ and the nature of the bound state in nuclear 



matter |jT8[ were also examined recently, but without a connection to free-space parameters. 
Here we focus on the use of an effective lagrangian with parameters fit from free-space and 
few-body data. 

The effective lagrangian consists of long-range interactions constrained by chiral sym- 
metry and the most general short-range interactions consistent with QCD symmetries. The 
coefficients of these short-range terms may eventually be derived from QCD, but at present 
must be fit by matching calculated and experimental observables in a momentum expansion. 
In free space, this expansion is in the external momentum over some breakdown scale A. 

For NN scattering, an EFT that consists only of short-range interactions is equivalent to 
the effective range expansion, which breaks down around A ca mT^/2. Using chiral symmetry 
to explicitly include pions in the lagrangian pushes this scale up to at least 300 MeV |T^-^, 



but it may be as high as the p(770) mass [p^-p^. This effective lagrangian can then be used 
to systematically predict other observables, including inelastic processes, up to the same 
scale. In the following, we discuss extending this procedure to nuclear matter. 

If the in-medium expansion parameter is simply the Fermi momentum kp over the free- 
space breakdown scale A, then A ^ 300 MeV implies that a useful EFT description of equi- 



libriuin nuclear matter (where kp ~ 260 MeV) might be impossible.^ The analysis of phe- 
nomenologically successful mean-field descriptions of nuclei from an EFT perspective reveals 



a systematic density expansion, with a characteristic mass scale around 600 MeV p5| , |26 
This suggests that a useful EFT expansion of nuclear matter is plausible if the free-space 
breakdown scale can be pushed as large as rrip, but a determination of how this expansion 
carries over to many-body systems is also required. 

Important questions about how to do power counting in a fully nonperturbative many- 
body calculation are as yet unresolved. Nevertheless, we can try to isolate some of the issues 
and techniques from free-space EFT analyses to explore how they translate to a many-body 
system. Our strategy here is to extend to the medium the perturbative matching example 
used by Lepage to illustrate EFT renormalization for ordinary two-body scattering [jl9| . 



This enables us to explore the convergence of non-scattering observables, the effective in- 
medium expansion parameter and breakdown scale, the use of error plots, and regularization 
dependence. Special considerations required due to the large scattering lengths characteristic 
of nuclear scattering are postponed to future studies. 

The convergence properties of an EFT in the medium can be cleanly tested by fitting 
to pseudo-data generated from an exact underlying potential that features a characteristic 
mass scale. A full EFT calculation must include long-range pion effects, but these do not 
modify the basic analysis presented here, so the role of the pion is ignored in order to keep 
the discussion as transparent as possible. Similarly, we omit higher partial waves, spin 
dependence, and relativistic corrections. For our discussion, the only key feature at finite 
density is the Pauli-blocking of occupied intermediate states. With a simple interaction 
between nucleons that is perturbative in the coupling, expansion parameters in free space 
and in the medium can be identified both analytically and using the corresponding error 
plots. 

In Sect. II we perform perturbative matching calculations in free space using a simple 
model potential, with both dimensional regularization and a cutoff regularization. We illus- 
trate the matching procedure in some detail, as it is likely to be unfamiliar to practitioners 
of conventional nuclear physics. In Sect. Ill we apply the EFT to finite density and show 
how the errors and breakdown scale extend to the medium. The implications for EFT 
calculations of nuclear matter are discussed in Sect. IV. 



II. MATCHING IN FREE SPACE 

Galilean invariance requires the interaction between two free nucleons of momentum ki 
and k2 to be independent of their center-of-mass momentum P = ki + k2. For clarity, we 
take the underlying potential to be a separable potential in the ^5*0 -channel^] that falls off 
at large momenta and depends on the relative momentum fc = |k| = i|ki — k2| only. 



^However, in practice numerical factors could improve (or degrade) the situation. 



^An extension to higher partial waves is straightforward |27| and does not introduce new features 
into the analysis. 



\^\V\.rVic\^l — T^jfio I 2\n 12 I 2\ ' \ ' 

M (k^ + 171'^ j[k'^ + m^) 

The mass m corresponds to the range (and non-locahty) of the potential; it plays the role of 
the underlying short-distance scale. We adopt an overall normalization oiAir/M, with M the 
nucleon mass. The dimensionless coupling a provides a perturbative expansion parameter. If 
a is 0{1), then the effective range parameters are of "natural" size, which means a constant 
of order one times, in this case, a power of m (see Appendix ^. 

An EFT can describe observables from the interaction in Eq. (|l|) to any desired accuracy 
as long as the details of the underlying potential are not probed. For scattering observables, 
this restricts the external momentum k to be much less than the characteristic mass m. The 
effective potential can then be written as a momentum expansion 

(k'ifEFTik) = Co + C,^-^^±^ + C y^p' + c S^'-p' + . . . , (2) 

which must be regulated. The two most common regularization schemes are dimensional 
regularization with power divergence subtraction (DR/PDS) |^ and cutoff regularization 
(CR) 0. In the latter case we use a gaussian separable cutoff in momentum space that 
strongly damps momenta above an arbitrary cutoff A^: 



{k'|l>EFT|k) 
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^_(fc2+fc,2)/2A2 _ j-g^i 



For perturbative calculations, analytic expansions can be obtained for these schemes in both 
free-space and in a uniform finite-density system. 

The idea of a perturbative matching calculation is that we can work to rfi^ order in the 
momentum expansion and determine the corresponding coefficients Co, ... , C2n order- by- 
order perturbatively in the coupling a by equating observables calculated from the "true" 
and EFT hamiltonians [|l^. We define the series expansion coefficients C2„ by 

oo 

C2n = E«'^2S- (4) 

s=l 

It is not necessary that perturbation theory for the observables converges in the low- 
momentum limit ,0 but for our example it does. We will carry out the matching for Co 
in some detail so that the generalization to higher orders and the extension to finite density 
is clear. We choose the on-shell T-matrix as our matching observable, since it has a natural 
perturbative expansion: the Born series T = V + VGqV + ■ ■ ■, with Go = {E — Hq)~^. 



^The main assumption is that the short-distance physics is perturbative. The coefficients will 
have convergent expansions in a since, as we wifi see explicitly, higher-order constants incorporate 
high momentum physics only and so are not sensitive to possible infrared divergences (e.g., from 
a long-range Coulomb potential). This means that perturbative matching is sufficient even if the 
constants will then be used in a nonperturbative calculation (see Ref. [^] for an example). 
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FIG. 1. Perturbative matching in free-space to 0{k'^). The double hne for the potential is a 
reminder that the potential is separable. 



We start our discussion with cutoff regularization since the physical content of the renor- 
malization is manifest. At 0{a), the on-shell T-matrix is just the on-shell potential V, 
which for k <^ {m, A^} can be expanded (it is sufficient to take k' = k for 5- waves): 
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where we have used the shorthand notation Oik"^) to denote natural corrections with the 
proper dimension multiplied by either k'^ /m? or /c^/A^ as applicable. Matching to this order 
fixes 
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and is indicated schematically in Fig. |l|a. 

Matching at O(a^) is illustrated in Fig. [pD. There are two contributions on the EFT 
side: (VGqV) evaluated using Cq fixed by Eq. (P), and {V) evaluated with the 0{a'^) 

(2) 

contribution to Cq, namely, Cq . The latter renormalizes Cq. To illustrate the physics of 
this renormalization, we consider explicitly the difference of (V^G'oV^)true and (V^GoV^)eft: 
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(7) 



Since we are only working to 0{k'^) in the effective potential, we can expand all k'^ depen- 
dence except for Gq, since q can be smaller than k: 
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The constant Cq is chosen to cancel the ^^-independent part of A(VG'oV), so that the net 
result for {V + VGoV) is 0{k'^): 

^° " Mm 2 .AFm ' ^^^ 



Choosing the cutoff Ac ~ m keeps the constants natural and gives the maximum range of 
validity for the EFT ^. 



(2) 



As expected from the uncertainty principle, the local interaction Cq is determined 



entirely by high-momentum (g > A^) contributions in the loop integral |T^ . This is because 
for q <^ Ac, the true and EFT tree-level results have already been matched in Eq. (^ to 
leading order in q^ and therefore cancel in the integrand of Eq. (P) . Note that Cq removes 
the Ac dependence while correcting for the contributions in the loop integral from the high- 
energy states. The ability to absorb the high-momentum components of the interaction 
into the constants of the effective potential is a generic feature of an EFT approach and is 
essential for a systematic prediction. 

If we add a long-range potential to the true theory, we must also reproduce its long- 
wavelength effects in the EFT, so we will still have agreement in the low-momentum region 
where q -C Ac, and the analysis goes through unchanged. If we extend the effective potential 
to include two constants, the 06*2 q^ piece serves to make the low- momentum part of the 
loop integrals agree to 0{q'^), and the high-momentum part of ^(yGoV) is absorbed into 
q;^(Co + C2 k"^)- As a result, the EFT reproduces the observables to 0{k'^). The addition 
of C2 requires an additional renormalization of Cq in cutoff regularization (denoted by 5Cq , 
see Appendix ^), which means that all the constants change with each successive order in 
the momentum expansion. 

Having convinced ourselves of the physical origin of the constants, we move to the simpler 
but physically less transparent dimensional regularization. It has been shown that power 
counting with only short-range potentials is regularization scheme independent for free-space 
observables p0|j28[] , as considered here. Comparing the expressions for the true and effective 
T-matrix to leading order in a again leads to 
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which fixes the first term in Eq. (H) to be 

^° -Mm' ^^^^ 

just as with the cutoff regulator. 

At the next order in a, the difference between {VGoV)tTue and (V^G'oV^)eft becomes 
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where the second integral is dimensionally regularized with D = 4 + e and in general is 
evaluated in the PDS scheme as |^ 
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Here /x is the DR renormalization scale. As before, the low-momentum parts of the in- 
tegrands already agree to 0{q'^), so the constant difference is from the high-momentum 
behavior. 

The results for the first two terms in the Born series are 
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which requires 
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for a proper match to 0{a'^,k'^). Note that for Ac = v^/^5 the DR/PDS constants are 
equivalent to the cutoff results Eqs. (H) and @. This agreement does not persist at higher 
orders in momentum. 

(s) 

Carrying the matching to one more order in a suggests a geometric series for the Cg 
in the DR/PDS scheme, which indeed sums to the full nonperturbative solution (see Ap- 
pendix 1^): 
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Extending the analysis to higher orders in momentum (or expanding the full results from 
Appendix ^ in powers of a) gives 
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As additional constants are added to the effective potential in DR/PDS, the previously 
fixed constants are not modified. This is not the case in CR (see Appendix 0) due to terms 
proportional to positive powers of the cutoff. 



The constant C4 requires additional input, since it does not contribute to the on-shell 
two-body T-matrix. It is tempting to match the true and EFT T-matrices off shell to 
determine C4, but this is never necessary, as only observables are required to fix the EFT 
constants. In particular, the additional on-shell constraint of matching the three-body 
scattering amplitudes can be used to find C4. Since we are in a regime where the 
coupling is perturbative, the Faddeev equations are simplified to a set of three-to-three 
tree-level amplitudes at second-order in the potential, 0{a'^). The result from matching is: 
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The second-order constant C4 as well as true three-body terms (needed to absorb the diver- 
gences in three-to-three loop- level amplitudes |^) do not enter nuclear matter calculations 
until higher order in the a expansion. 

By construction, the EFT systematically reduces the error order-by-order in the mo- 
mentum expansion. We can see this explicitly by examining the truncation error in (T) to 
(9(a;^), which is just {V + VGqV) , between calculations using the true and EFT potentials.^ 
For an effective potential only containing Co, the error is (9(fc^) 
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and adding C2 brings the error to 0{k'^) 
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Note that the truncation error with CR depends on Ac while in the DR/PDS scheme the 
truncation error is independent of //. Extending the analysis to nonperturbative calculations 
and finite density will in general require numerical solutions, and so a connection between 
the analytical results above and a graphical error analysis is important. The error plots 



introduced by Lepage [19| are useful in this regard 



^For the present perturbative case, this is more convenient than looking at the error in kcotS, 
which is appropriate for a nonperturbative calculation. In the error plots, we consider only the 
real part of the difference in {V + VGqV). The imaginary part follows from unitarity, which the 
EFT reproduces order-by-order. 
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FIG. 2. |ARe(k|f|k)| to 0{a^) vs. k/m for a = -1/2 and with M = 940 MeV and 
m = 600 MeV. Results for cutoff regularization (CR) with Ac = m are shown on the left and 
for dimensional regularization (DR/PDS) on the right. Note that Ug = — 2/3m and r^ = 9/m. 



Keeping constants to 0{k^), 0{k^), and 0{k'^) in the effective potential Eq. (Q) or Eq. (^ 
leads to successively better approximations to {V + VGqV), as seen by the analytic expres- 
sions for the error Eqs. (|20D and (pT]) and shown graphically in Fig. ^ With each additional 
order, the slope of the error increases, reflecting the improved truncation error. (If a long- 
range potential is added to both the true and effective potentials, the absolute error in the 
DR/PDS power-counting scheme will decrease but the error is always 0{k'^) |2yj2^.) For 
k <^ m the error is dominated by the first unmatched term, so the lines are straight on a log- 
log plot. The EFT should break down when the external momentum probes the details of 
the underlying potential; graphically the intersection of the line^ indicates the approximate 
breakdown scale A. In our example, we see that A ^ m, as expected for a natural theory 
[with a ~ C)(l)]; indeed, all the errors are of the same order for k ^ m. This breakdown 
scale does not change as more orders in a are included in the contact interactions, but the 
accuracy does improve. 



III. FINITE DENSITY OBSERVABLES 

Now that the constants are determined to some fixed order in the momentum and cou- 
pling expansion, the EFT can be used in a many-body calculation and compared with 
results for the true potential. Since the EFT constants are matched to data in free space. 



^To determine the intersection scale, one should extend the lines from the straight regions {k <^ m) 
rather than looking at the actual intersection. 




FIG. 3. Integration regions F and F. 



renormalization-schenie independence requires a comparable level of sophistication for a 
calculation in the medium. If the free-space T-matrix is solved using the nonperturbative 
Lippmann-Schwinger equation, the in-medium T-matrix (denoted by F) will at least require 
a nonperturbative solution of the Bethe-Goldstone equation ||30|. However, this leads to 
many complications at finite density. For example, the two-nucleon Green's function G now 
depends on the interaction between the nucleons, which itself depends on F, leading to a 
self-consistent equation. The integrations over nucleon states are much more complicated 
since occupied intermediate states are Pauli-blocked. Furthermore, additional terms may 
also be required, such as interactions between particles and holes. The development of a 
consistent yet tractable power counting scheme for finite density is a critical issue for future 
investigations. 

Here, however, these issues are greatly simplified by virtue of having perturbative inter- 
actions and working only to second order in a. Replacing G with the free-space propagator 
Go only changes F at 0{a^), and the perturbative calculation including the Pauli exclusion 
projection operator can be performed analytically to this order as well. The equation for F 
can be written in momentum space as 
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with the region of integration F given by the exclusion of the two Fermi spheres depicted in 
Fig. |: 

Jr 27r 3 J 27r 3 V 2 ^ ^/ V 2 ^ 
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The subscript on Fp in Eq. (|2^ ) is a reminder that, due to the presence of the Fermi sea, 
there is an explicit dependence on the center-of-mass momentum P of the two interacting 
nucleons, which to this order appears only in the integration limits. 

The matrix elements of Fp can be used to find the ground state energy per particle for 
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FIG. 4. Perturbative matching for the energy in nuclear matter at 0{k^p) 
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Here g = {2S + 1)(2/+ 1) is the spin-isospin degeneracy factor, p = gkp/ {6n'^) is the density 
of nucleons, and the integration region F is the intersection of the two Fermi spheres, as 
depicted in Fig. |^. We will use Eg,s./A as our basic finite density observable. 

We start with the DR/PDS effective field theory. Evaluating Eq. (|^) to 0{a) is just 
the Hartree-Fock approximation; the relevant diagrams are shown in Fig. ^. Both the EFT 
and true results can be analytically integrated and expanded (see Appendix y) to give (with 
^ = kF/m) 
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with terms only up to C4, shown in the EFT calculation. As each new constant C2n is added 
to the effective potential and the values found in free space are used from Eqs. (pJ| ), (|15|), 



(p!7|), and (p!8D, the EFT result for i?g.s./A coincides with the true result up to truncation 
errors of 0{k'^p^^) and 0{a'^). There is no surprise here: The perturbative matching at finite 
density must work to this order since the EFT is merely parametrizing the Taylor expansion 
of V^rue before the integration over the intersection of the Fermi spheres in Eq. (pS)). 

Taking the analysis to the next order in a requires an evaluation of the Pauli-blocked 
{VGqV) shown in Fig. P^, and hence is a better test of how the EFT carries over to nuclear 
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matter. The low-momentum part of the true and EFT integrals still agree up to (9(/c^"), 
since the non-Pauli-blocked part includes the same momentum states in each. If kp <C m, 
the Pauli blocking has no effect on the high-momentum parts of the integrals, which are 
therefore the same as in free-space. Thus the renormalization in free space of C2J by C2J 
carries over directly. Integrating over F to find the energy leaves the same truncation error, 
0{Pp^^), as at order a. 

The steps involved in actually carrying out the integrations are outlined in Appendix ^ 
and the final expressions are given there. If we expand the energy per particle in powers of 
^, the true result is 
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and the EFT result using DR/PDS regularization is 
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Again, substitution of the expressions up to C2n from Eqs. (|n]), (|T3|), and (|T7|-|T^) removes 
the /i dependence and leads to agreement with the true result up to truncation errors of 
0{k'^p^^). For a perturbative DR/PDS calculation with kp <^ rn^ these errors are simply 
given by the first unmatched term in the expansion of the true result. For one constant the 
leading error is 
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The full error is depicted graphically in Fig. |^, which shows a similar breakdown scale to 
that in free space. 
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FIG. 5. Error in the ground state energy per particle in nuclear matter to 0{a'^) vs. ^ = kp/m 
for a = —1/2 with M = 940 MeV and m = 600 MeV. The error is divided by S,^ for clarity. Results 
for cutoff regularization (CR) with A.c = m are shown on the left and for dimensional regularization 
(DR/PDS) on the right. 

We emphasize that the contribution from the term in Eq. (pSD proportional to Cl is 
needed to reproduce the ^^ term in Eq. (P7|). Since this coefficient does not contribute to 
on-shell two-body scattering, it might be interpreted as an "off-shell ambiguity." However, 
as discussed above, it is determined by on-shell three-body scattering, so there is no need 
to consider off-shell effects in the two-body sector. At higher order in a, we will need full 
three-body contact terms, matched in the three-body sector, to account for high-energy 
contributions [^]. 

We find analogous results with cutoff regularization. The expanded 0{a) energy is 
(keeping terms to 0{k'^) in the effective potential): 



CR: '"* 
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(9 



am 



(Mm 



ci^'^e + ^m^ 






e 



+ 0{C) (31) 



(2) 



and using the explicit expression for 5Cq from Appendix |B|, the Oi^a^) result is 
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(32) 



These expressions parallel those found in the DR/PDS scheme in Eqs. ( p6D and (0). Tak- 
ing the values for the constants found from matching in free-space (see Appendix 0), and 
substituting them into Eqs. (^) and (^) produces Ac-independent results at each order in 
the momentum expansion (although the truncation error depends on A^). The leading error 
for one constant is 
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e + o{e) , (33) 



and for two constants it is 
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(34) 



Just as in the DR/PDS case, the error decreases by a power of ^^ for each additional order, 
as shown in Fig. |^. This again shows the breakdown scale at finite density is similar to the 
breakdown scale in free space. 

In summary, at finite density the energy per particle has similar error plots to those in 
free-space when AE^^^JA is plotted against the Fermi momentum kp- The improvement is 
by two powers of kp with each additional order (although both even and odd powers of kp 
are present in the energy per particle). The apparent breakdown scale is still K k. m, without 
any numerical factors significantly different from unity. The important observation is not 
the precise breakdown scale, which depends on the details of the underlying interaction, but 
that the breakdown in free space and in the medium are closely related. The analysis is 
easily extended to higher partial waves for both regularization schemes, which verifies the 
matching and the conclusions about the breakdown scale. 

IV. DISCUSSION 



The model we have studied shows that an EFT expansion for NN scattering leads to 
a low-density expansion of the energy per nucleon in nuclear matter. In particular, in the 
former one finds scattering observables order-by-order in the nucleon relative momentum 
over a breakdown scale . In the latter, one finds the energy per particle order- by-order in 
the Fermi momentum kp over essentially the same scale. We have checked that this basic 
conclusion is independent of the details of the model potential. We expect that this type of 
expansion will also result when long-distance pion physics is explicitly included, in analogy 
to the modified effective range expansion in free space ||20[| . 

Of course, since this is a perturbative calculation it can only be suggestive. We will need 
to clarify a variety of interrelated issues: nonperturbative power counting in the medium, the 
impact of pions, the elimination of off-shell ambiguities, the role of three-nucleon interactions. 
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However, if we assume the basic results here carry over to a full nonperturbative treatment, 
we can speculate on how well a leading order (LO) and next-to-leading order (NLO) EFT 
calculation will do in reproducing nuclear matter saturation. 

In particular, we can use phenomenological energy functionals of nuclear matter (that is, 
the energy expressed as a function of the density) to anticipate what a full EFT calculation 
might achieve. Nonrelativistic and relativistic mean-field functionals provide parameteriza- 
tions of the energy per particle near the equilibrium density, with individual contributions 
characterized by a density expansion parameter of about 1/5 [^. These functionals are fit 
directly to finite density data. As such, they include (approximately) the effects of three- 
body forces and short-range correlations, even though the formalism is that of Hartree or 
Hartree-Fock |3^. (There are certainly important limitations in the descriptions; for exam- 
ple, there is no kjp term in nonrelativistic mean-field energy functionals.) 

If we assume that the expansions about kp = are reasonably determined by fits near 
(^f) equilibrium < A, then wc cau usc them to predict how far in a truncated kp expansion we 
need to go to see saturation]^ Ordinary Skyrme Hartree-Fock functionals are truncated at 
p^, which is p^ oc kp in the energy per particle. As more terms are added to the functional, 
it should become a more reliable "extrapolation" of an expansion about kp = 0. So we 
can look at a generalized Skyrme model, fit out to kp but then truncated at kp, to see 
how saturation changes. In this case, the equilibrium density decreases by about 20% while 
the binding energy decreases by roughly 5 MeV [Q. There is an energy minimum when 
only terms up through kp is kept, but it is far removed from the empirical point. Thus we 
expect a qualitative reproduction of saturation if we keep terms up through kp, as long as 
the effective A is sufficiently high. 

A similar conclusion is reached by examining relativistic mean-field point-coupling mod- 
els. These models contain contributions to all orders in the kp expansion. One finds that 
a saturation minimum is first achieved with an expansion truncated at kp but that the kp 
term is needed to get reasonably close. The quantitative results again suggest that an EFT 
calculation up to kp might reproduce nuclear saturation with the equilibrium binding energy 
given to within about 5 MeV and the density to within about 20%. 

The mass scale that sets the convergence of contributions to mean-field energy functionals 
is of order 600 MeV ||32|. However, this is not necessarily the same as the breakdown scale. 
If one constructs error plots for relativistic mean-field functionals by comparing the full 
energy as a function of kp to the same functional truncated at kp, at kp, and so on, a 
breakdown scale for the expansion about kp = is identified. For point-coupling models, 
this breakdown scale is close to 400 MeV but for models with meson fields it is actually 
below the equilibrium kp [Q. This discrepancy is not yet understood. 

One might ask how saturation is realized in such a density expansion. The existence 
of a minimum in the energy per particle suggests that different orders in the expansion 
must be playing off each other; a natural conclusion is that this should only occur near the 
breakdown scale (where all terms contribute equally) and not at kp <^ A. The mean-field 
phenomenology suggests a resolution: there is an interplay of different orders but it is highly 



^A nonperturbative EFT calculation would include contributions to all orders in kp and therefore 
may do even better. 
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restricted. In particular, repulsion from the kp and kp terms becomes important compared 
to the attraction from the kp piece well below the breakdown scale A. Furthermore, this 
happens because the coefficient of the kp term is unnaturally small;[] in relativistic models 
this is a direct result of cancellations between Lorentz scalar and vector contributions that 
are each of natural size. The cancellations leading to a small kp term do not recur in higher 
orders, so further contributions just correct the equilibrium properties and do not change 
the qualitative behavior. 

The convergence of the density expansion is also consistent with the convergence of 
a conventional expansion of the Galitskii equation in powers of kp times the radius of a 



strong, short-range (hard-core) potential |^, which we can associate with 1/A. In the 



effective field theory treatment, this short-distance expansion parameter arises naturally 
when long-distance pion physics is treated explicitly. In free space, consistently including 
the pion takes us from an effective range expansion with a breakdown scale of order the 
pion mass to a modified effective range expansion with a breakdown scale set by the next 
important mass scale. We expect an analogous situation at finite density. 

So what are the ingredients of a minimal EFT calculation of nuclear matter saturation? 
First, we need the breakdown scale for NN scattering in free space to be as large as possi- 
ble. There are indications that a complete NLO calculation using cutoff regularization and 
including two-pion irreducible contributions (and possibly the NNLO two-pion pieces) may 
achieve A ^ 600 MeV ||23| . Then we must include all contributions that will generate terms 
in the energy per particle up to kp. This implies including two constants in each S-wave, 
one in each P-wave, and an S-D mixing term. One- and two-pion exchange should be in- 
cluded in all waves. Finally, a three-body contact term will be needed ^. Note that this 
recipe includes the standard ingredients of nuclear matter phenomenology: the pion tensor 
force, velocity dependence of the interaction, mid-range attraction from two-pion exchange,^ 
and three-body forces. While much remains to be done before a consistent nonperturbative 
calculation is available, we are optimistic that our speculations will be tested in the near 
future. 
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^Here "small" means roughly half the size one might expect; it is not fine-tuned to be close to 
zero. If the coefficient is doubled in size, equilibrium does occur at fci? ps A and bound by 200 to 
300 MeV. 

^Note that the NLO or NNLO two-pion potential does not include contributions from two pions 
interacting "in flight," which is counter to the usual intuition from nuclear phenomenology about 
the mid-range attraction. 
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APPENDIX A: NONPERTURBATIVE FREE-SPACE SOLUTIONS 

Scattering from the S-wave separable potential 

(k iV^truelk) = -Jf ^k^^m^)ik- + m^) - ^o(^' ^ ) ' (Al) 

can be solved analytically to give 

T,-\k, k) = V,-\k, k) + ^ l^tk - ^^^ (A2) 

or 

Air k'^ — m^ , , , 

Expanding k cot (5o in powers of k"^ gives the effective range expansion for this potential: 

A:cot5o = -m(\ + -) + -(\--]e- ^fc^ . (A4) 

\1 a) m \2 a/ am-^ 

The effective-range parameters are natural (i.e., characterized by the underlying mass scale) 
for a = C(l). 

Regularizing the EFT with DR/PDS gives 

which can readily be expanded in momentum and solved for the S'-wave low-energy constants 

C2n'- 



c, 



a ZTcm \ 4/ Cq a Attitl-^ 



2 m 



Mm 
Expanding these results in powers of a gives the results in Eqs. (0), (P^, and ([TtHISD. 



Matching the tree-level on-shell three-body scattering amplitude IQ to leading order in a 
gives C*4 = C^jra^. 

APPENDIX B: TWO CONSTANTS WITH CUTOFF REGULARIZATION 

For two terms in the effective potential, the expression for Cg [Eq. (|^)] gets an extra 
piece: 

K'? = ^5^ (Cfcl" + 5(C»Va?) , (Bl) 

with 

The C(a^) piece of C2 has more terms compared to the DR/PDS result Eq. ([T7|): 
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APPENDIX C: NUCLEAR MATTER INTEGRALS 



For an even-parity potential, Eqs. ( [2^ ) and ( [25| ) give the first-order contribution to the 
energy per particle as 



E^ll 9(9-1) f d'Pd'k ,^ ,^^^,^ ,_9{9-h 4 



2p Jf {2tt 



{k\V\k) 



Ti^p 



m X . 



(CI) 



For an odd-parity potential, the degeneracy pre-factor becomes g{g + 1). Introducing scaled 
momenta s = ^Y* /kp and t = k/fc^, then 



l = i^ I dss^ I dx I dtt^J^it) 



JO 



with 



Z±(X) = ±SX + \/l — 5^(1 — X 



T{t) = kl {kFt\V\kFt) 



(C2) 



(C3) 



Integrating by parts with respect to x and then s, and changing variables simplifies the 
integral to 



6 7o 
which for the potential in Eq. (|l|) gives 

27rm 
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(C4) 



(C5) 



The integral can also be directly evaluated for the DR/PDS and CR effective potentials. 
The second-order contribution El^l/A follows from the same analysis as above, but with 

J^(t) = klV -——ikFt\V\kFVL)- -{kFn\V\kFt) = dy duu'^n{t,u) , (C6) 



where 
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We divide the integrals into two parts: 

X = Xi — X2 , 
with 

rl rl pz—ix) rOQ 
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(C9) 
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The first integral Xi can be evaluated as in Eq. ( |C4|) . For the potential Eq. (|l|) , this leads to 
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The second integral X2 is not conveniently expressed in closed form. For the potential Eq. (^, 
it can be expanded in ^ to give 
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f m 



V105M 

C^^ /26101 
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Results for the EFT are found similarly. 
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